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Abstract 

Based on a three-potential formalism we propose mathematically well- 
behaved Faddeev-type integral equations for the atomic three-body prob- 
lem and descibe their solutions in Coulomb-Sturmian space representation. 
Although the system contains only long-range Coulomb interactions these 
^ equations allow us to reach solution by approximating only some auxiliary 

• ■ short-range type potentials. We outline the method for bound states and 

O ' demonstrate its power in benchmark calculations. We can report a fast con- 

c/2 , vergence in angular momentum channels. 

P^l The Faddeev equations are the fundamental equations of the three-body problems. Be- 

sides giving a unified formulation, they are superior to the Schrodinger equation in many 
^ ■ respects: in incorporating the boundary conditions, in treating the symmetries, in handling 

the correlations, etc.. Nevertheless, their use in atomic three-body calculations is rather 



O i scarce [0-^, and these calculations showed that to reach a reasonable accuracy many chan- 

r^ I nels are needed. So, the behef spread in the community that the Faddeev equations, the 

O ' fundamental equations of three-body systems, are not well-suited for treating atomic three- 



body problems and other techniques can do a much better job, at least for bound states. 
In Refs. Ill-Q the Faddeev equations were used in such a form that, the solution could 

O i be reached only by some kind of approximation on the long-range Coulomb potential, and 

c/3 I thus the convergence were ensured only via the square integrability of the bound-state wave 

r^- function. 

Ph' The aim of this paper is to solve the atomic three-body problems by approximating only 

short-range type interactions. We invoke a newly established "three-potential" formalism 
k> ! and derive such kind of Faddeev-type integral equations which contain only short-range type 

5^ I interactions as source terms. We solve the equations by approximating only the short-range 

type source terms. We point out that although we are working with finite matrices the 
wave functions possesses correct three-body Coulomb asymptotics. Finally, as compulsory 
benchmark cases, we calculate the helium atom, the positronium ion and the muonic hydro- 
gen molecule ion, and will observe a fast convergence with respect to angular momentum 
channels. 

The "three-potential" formalism was designed for solving nuclear three-body problems 
in the presence of Coulomb interaction. The method was presented first in bound-sate 
calculations |^ and was extended to below-breakup scattering calculations P where the 

1 



notion of the "three-potential" formahsm was also introduced. In this solution to the nuclear 
three-body problem, in the spirit of the two-potential formalism, all long-range interactions, 
i.e. all Coulomb interactions, are put, a la Noble 0, into the "free" Green's operator. 
Invoking again the two-potential formalism, the incalculable Faddeev- Noble Green's operator 
was linked to a simpler channel distorted Green's operator. The "three-potential" formalism, 
thus consists of repeated applications of the two-potential formalism, and results in integral 
equations which contain only short-range interactions as source terms. 
The Hamiltonian of an atomic three-body system reads 

H = H' + v^ + v'^ + v^, (1) 

where H^ is the three-body kinetic energy operator and v^ denotes the Coulomb interaction 
in subsystem a. We introduce here the usual configuration-space Jacobi coordinates ^a and 
fja] ^a is the vector connecting the pair (/5,7) and fja is the vector connecting the center of 
mass of the pair (/?,7) and the particle a. Thus v^ only depends on ^q,, i.e. v^ = f^ (^a)- 

The peculiarity of Hamiltonian (|l|) is that all potentials are of long-range type, and 
thus the procedure of Refs. [^Q cannot be applied without any modification. However, 
the physical role of the Coulomb potential is twofold. Its long-distance part modifies the 
asymptotic motion, while its short-range part strongly correlates the two-body subsystems. 
So, we may split the Coulomb potential as a sum of long-range and short-range terms, 

i;^=t;W+t;W. (2) 

Short-range and long-range interactions play entirely different roles in integral equations. 
While short-range interactions contribute to the source terms, long-range interactions should 
always be put into the Green's operator. Following Faddeev's procedure §], we split the 
wave function into three components, 

\^) = m + \i^p) + \i^,), (3) 

and, for the components, we arrive at the Faddeev-Noble-type equations 

lV'a) = G«(E)[.;W|^,)+^i^)|^,)], (4) 

with a cyclic permutation for a,/5, 7. The Green's operator is defined as 

G« (.) = (. -/7° -.«-.«-<)-. i^))- 

= {z - H' - v^ - v^^ - v\'^r\ (5) 

In the spirit of the "three-potential" formalism, we relate this Green's operator to the simpler 
channel distorted Green's operator via the resolvent relation 

(6) 



(7) 







G^\z) = 


--G^{z) + G^{z)U"G'f>{z\ 


with Ga{z) 


andf/" 


defined by 








Ga{z) ■■ 


= {z-H^-v^^-u^l^r 



and 

t/" = 4' + <')-«i^ (8) 

respectively. Here we have introduced the auxihary potential m^'-* which is acting in coordi- 
nate T]a and is required to have the asymptotic form 

Ua\Va) ~ e„(e^ + e^)/?7„ as ?7„ ^ cx). (9) 

In fact, it is an effective long-range interaction between the center of mass of the subsystem 
a (with charge e/3 -|- e^) and the third particle (with charge Cq,). The importance of u^^ lies 
in the fact that it asymptotically compensates the Coulomb-tail of u^ + f^'-' in Eq. (H), and 
thus U" decays faster than the Coulomb interaction in the two-body channel a. It does 
not influence the character of the asymptotic motion, so, although it may contain terms 
which behave like l/rja as rja — > cxd, it can be treated as short-range interaction. So, in Eq. 
(I) and in Eq. (H), which constitute the basic integral equations of the atomic three-body 
problem mediated by the "three-potential" formalism, only short-range type interactions 
are appearing as source terms, thus the equations are mathematically well behaved. They, 
similarly to the nuclear Coulomb three-body problem [§,|1, can be solved in Coulomb- 
Sturmian (CS) space representation. 

The CS functions, which are the Sturm-Liouville solutions of the hydrogenic problem 0, 
are defined in configuration-space for some angular momentum state / as 



{r\nl) = \J ^^^2l + l)\ i^bry'''^-'^L':^\2br), (10) 

where n = 0, 1, 2, . . ., L represents the Laguerre polynomials and 6 is a fixed real parameter. 
With the functions {r\nl) = {r\nl)/r they form a biorthonormal basis. 

Since the three-body Hilbert space is a direct sum of two-body Hilbert spaces, the ap- 
propriate basis in angular momentum representation may be defined as a the direct product 

\niylX)a = [\nl)a(S)\iyX)a]L, (n,z/ = 0,1,2,...), (11) 

with the CS states of Eq. (|T0|). Here / and A denote the angular momenta of coordinates ^ 
and 7], respectively, and they are coupled to the total angular momentum L. Now the com- 
pleteness relation takes the form (with angular momentum summation implicitly included) 

TV 

1 = lim y^ \niylX)a a(nz//A| = hm 1^, (12) 

where {^aVa\^^^^)a = ^/C,aVa ('Ca^al nvlX) a- It should be noted that in the three-particle 
Hilbert space we can introduce three equivalent bases belonging to fragmentation a, (3 and 

7- 

In equations (^) we make the following approximation: 

|^„) = G«(E)[1>W4|^^) + irfl^lV^,)], (13) 



i.e. we approximate the short-range potential vj^^ in the three-body Hilbert space by a 
separable form 

N 

vi'^- E \niylX)av!alp{n'u'l'X'\, (14) 

n,i/,n' ,u'=0 

where i^*j = (1 — Sap) Q,(ni//A|f^'^^|nV/'A')/3. In Eq. (|l^ the ket and bra states are defined 
for different fragmentations depending on the environments of the potential operators in the 
equations. 

Multiplied by the CS states ainiy^M from the left, Eqs. (|T3[) turn into a linear system of 
homogeneous equations for the coefficients of the Faddeev components ip = Q,(rai/ZA|^a): 

{[G!^'\E)]-'-v^''^}^ = 0, (15) 

with 

G« = „(m4A|G«|n'iW^A%. (16) 

Eq. ( |15|) is solvable if and only if 

det{[G'-'\E)]-' - v!^'^} = 0. (17) 

Notice that the matrix elements of the Green's operator, which contains all long-range 
terms, are needed only between the same partition a whereas the matrix elements of the 
short range potentials occur only between different partitions a and (3. The latter can be 
evaluated numerically by making use of the transformation of Jacobi coordinates |TU|. 



The matrix elements (|I6|) can be obtained by solving the Eq. (|g), which is, in fact, a 
two-body multichannel Lippmann-Schwinger equation. If we perform again the separable 



approximation (IT^) on potential f/°, with the help of the formal solution of Eq. (Rf) we may 



now express the inverse matrix (Gl/{E)) ^ as 

(G«)-^ = (GJ-i-I/", (18) 

where 

G„ = o.{ru4X\Ga\n'^'X')a (19) 

and ir = o.{nvlX\m\n'v'l'X')^. 

While the latter matrix elements may again be evaluated numerically, for the calculation 
of the matrix elements in Eq. (|1^) we proceed as follows. Since we can write the three- 
particle free Hamiltonian as a sum of two-particle free Hamiltonians, H^ = h2 + /i2 , the 

channel distorted Green's operator Ga of Eq. (|^) appears as a resolvent of the sum of two 
commuting Hamiltonians h? = h2 +v^ and /il'^ = h'l +uii\ which act in different Hilbert 
spaces. Thus, according to the convolution theorem |jlT|] we can express the three-body 
Green's operator Ga by an integral of two-body Green's operators 

Ga{z) = {z- hi - h^ir = ^ £ dz' {z' - hl)-\z -z'- h^-\ (20) 



Here the contour C should go, in counterclockwise direction, around the continuous and 
discrete spectrum of W^ in such a way that the resolvent of h^S^ is analytic in the domain 
encircled by C [cf. Fig. |I]]. This condition can only be fulfilled if we define the auxiliary 
potential u'^^ in such a way that the discrete spectrum of /i^2 does not penetrate into C. 
The matrix elements (IT9|) can be cast into the form 



G{E) = ^.idz' M\{z' - hl)-^\n'V)^ ^{u\\{E - z' - /i«)->'A%, (21) 

where both matrix elements occurring in the integrand are known from the two-particle case 



0- 



After solving Eq. ([TsD for the coefficients ^l) we can construct the Faddeev components. 
Considering Eqs. (|13|), (^ and (^) we get: 



1^.) 



^ £rf/ (.' - hlYM)^ [E-z'- hfx>^). 



e„, (22) 



where C_^ = {G^) ^{G:^') ^[v^a^ipg + v^^ip ], and the functions in the integrand are also 
known from two-particle case [|T^]. So, in this formalism, the Faddeev components appears 
as a linear combination of a convolution integral of Coulomb-like functions. Since in this 
procedure we have approximated only short-range type interactions {ipa) should possess 
correct three-body Coulomb asymptotics. 

To demonstrate the efficiency of the method we present a few benchmark calculations. 
As cut-off functions we use the exceptionally smooth error functions and define the auxiliary 
potentials as: 

v^\U=^^K^Uv^iU, (23) 

vl'\U=^M^UvaiU (24) 

and 

u^lKVa) = Ae"'^''' + e„(e/3 + e^) erf(cjr/«)/r/«, (25) 

where u, A and k are parameters. In defining u^^\ in order to prevent the penetration of 
the bound-states of /i^2 i^^o ^^e contour C, we have added a repulsive gaussian term. 

We examine the convergence of the three-body energy with increasing N, the number 
of terms in the expansion of the short-range potentials v^^ and f/". In Table |I| we present 
results for the binding energy for helium atom (with infinitely massive core), positronium 
ion {e~e~^e~) and the muonic hydrogen molecule ion {ppfi~)- For comparison we give the 
corresponding result of Faddeev calculations of Ref. 0,[13|, where channel-by-channel com- 
parison is possibble, and of Ref. [Q , where this kind of comparison is not possibble, and also 



of results of a very accurate variational calculation [T^\. In all cases we can observe that, 
similarly to what we have experienced with nuclear potentials without and with Coulomb 
interaction ^ , convergence up to 6-7 significant digits is comfortably achieved with terms up 
to A^ ~ 19 applied for n and z/, and the results are in good agreements with well established 



benchmark results. We can see that with an appropriate choice of cj, a rapid convergence 
is reached in the partial wave expansion. The convergence in angular momentum channels 
is much faster then in earlier works 0]. To overcome the poor convergence, the solution of 
the Faddeev equations in total angular momentum representation was proposed [0,|I, which, 
however, results in three-dimensional equations. In the light of our results, this seems to be 
superfluous, since the poor convergence is a consequence of the ill-behavior of the applied 
form of the Faddeev equations. 

It should be noted that in the solution of the Coulomb three-body problem the splitting 
of the Coulomb interaction into short- and long-range terms were already suggested before 



r5| . However, this splitting were performed in the three-body configuration space and the 
resulted Faddeev differential equations were used only for scattering state calculations. In 
our case the splitting is performed in the subsystem (two-body) Hilbert space and can be 
effective also for bound-state problems. 

We have proposed a set of Faddeev and Lippmann-Schwinger equations for atomic three- 
body systems using the newly established "three-potential" formalism. The Coulomb in- 
teractions were split into long-range and short-range terms and the Faddeev procedure was 
applied only to the short-range potentials. The resulting modified Faddeev equations are 
mathematically well behaved since the source terms are of short-range type and all the long- 
range interactions are kept in the Green's operator. This Green's operator was calculated 
from the channel distorted Green's operator via two-body multichannel Lippmann-Schwinger 
equations whose kernels are also of short-range type. Because the channel distorted Green's 
operator appears as a resolvent of the sum of two commuting two-body Hamiltonians, it can 
be represented as a convolution integral of the corresponding two-body Green's operators. 
The use of Coulomb-Sturmian functions is essential as it allows an analytic representation of 
the two-body Coulomb Green's operator on the complex energy plane, which makes straight- 
forward the calculation of the convolution integral and the incorporation of the infinitely 
many bound states in an attractive Coulomb potential. In this way the atomic three-body 
problem, which contains only long-range Coulomb interactions, were solved through approx- 
imating only some auxiliary short-range potentials. From this follows that the wave function 
possesses correct asymptotics. The extension of this formalism to below-breakup scatter- 
ing calculations is analogous to the nuclear Coulomb case 0. We hope, that these unique 
advantages of this method allow us to extend its scope to above-breakup calculations. 

This work has been supported by OTKA under Contracts No. T17298 and T020409. 
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TABLES 



TABLE L Convergence of the binding energy of the He atom, the positroniuni ion (e~e~^e~) 
and the niuonic hydrogen molecule ion ppfj.~ with increasing basis for the separable expansion 
taking into account angular momentum states up to / = 0, / = 2, / = 4, / = 6 and I = 8. N 
denotes the maximum number of basis states employed for n and u. The values for He atom and 
positronium are given in atomic units {h = rUe = e^ = 1) and the values for ppn~ molecule ion are 



given in modified atomic units (h = m^ = e = 1). 






Angular momentum channels 
N /=0 1=2 1=4 1=6 


1-- 


= 8 


He atom: 6 = 5, w = 1.5, A = 5 and k = 2 



2.892973 


2.903778 


2.903733 


2.903724 


2.903723 


2.892976 


2.903780 


2.903735 


2.903726 


2.903725 


2.892976 


2.903781 


2.903736 


2.903727 


2.903725 


2.892976 


2.903781 


2.903736 


2.903727 


2.903726 


2.892977 


2.903781 


2.903736 


2.903727 


2.903725 


2.892977 


2.903781 


2.903736 


2.903727 


2.903725 


2.892977 


2.903781 


2.903736 


2.903727 


2.903725 


2.892977 


2.903781 


2.903736 


2.903727 


2.903725 



Faddeev in Ref. @ 
with L up to 20 : 




2.9076746287 2.9044691278 


- 


2.9037959292 
2.9037224647 


Variational in Ref. H 








2.903724376984 




e 


'e"*"e~ ion: b = 0.8, u) = 0.2, A = 5 and n 


= 0.3 





0.250602 


0.261753 


0.261985 


0.261997 


0.261998 


0.250629 


0.261761 


0.261992 


0.262005 


0.262006 


0.250621 


0.261763 


0.261994 


0.262006 


0.262008 


0.250572 


0.261758 


0.261990 


0.262002 


0.262003 


0.250595 


0.261762 


0.261993 


0.262005 


0.262006 


0.250562 


0.261760 


0.261991 


0.262003 


0.262004 


0.250569 


0.261760 


0.261991 


0.262004 


0.262005 


0.250559 


0.261760 


0.261991 


0.262003 


0.262005 



Faddeev in Ref. || 














0.2620217 


Variational in Ref. [p^ 














0.2620050702314 




pp^i 


ion: b = 


= 2.5, u = 


= 0.8, A = 


= 20 and k = 


= 0.8 





12 
13 
14 
15 
16 
17 
18 
19 



0.462815 
0.458168 
0.455523 
0.453833 
0.452783 
0.452065 
0.451559 
0.451202 



0.486698 
0.486700 
0.486700 
0.486695 
0.486695 
0.486695 
0.486695 
0.486695 



0.492922 
0.492921 
0.492922 
0.492918 
0.492918 
0.492918 
0.492918 
0.492918 



0.494063 
0.494063 
0.494063 
0.494060 
0.494060 
0.494060 
0.494060 
0.494060 



0.494309 
0.494309 
0.494309 
0.494306 
0.494306 
0.494305 
0.494305 
0.494305 



Faddeev in Ref. ||4| 



0.4943867 



Variational in Ref. m4 



0.4943867 
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FIG. 1. Contour C for the integral for Ga{E) in case of the three-body bound-state problem 
with attractive Coulomb interaction. The contour C encircles the continuous and discrete spectrum 
of h9 and avoids the continuous and discrete spectrum of Wql . The discrete spectrum of h9 and 

hr^l is denoted by dots and crosses, respectively. 



